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ABSTRACT 

We investigate the photometric properties of the early type Brightest Cluster Galaxies 
(BCGs) using a carefully selected sample of 85 BCGs from the C4 cluster catalogue 
with redshift less than 0.1. We perform accurate background subtractions, and surface 
photometry for these BCGs to 25 magarcsec -2 in the Sloan r-band. By quantitatively 
analysing the gradient of the Petrosian profiles of BCGs, we find that a large fraction 
of BCGs have extended stellar envelopes in their outskirts; more luminous BCGs 
tend to have more extended stellar halos that are likely connected with mergers. A 
comparison sample of elliptical galaxies was chosen with similar apparent magnitude 
and redshift ranges, for which the same photometric analysis procedure is applied. We 
find that BCGs have steeper size-luminosity (R oc L a ) and Faber- Jackson (L oc 
relations than the bulk of early type galaxies. Furthermore, the power-law indices (a 
and (3) in these relations increase as the isophotal limits become deeper. For isophotal 
limits from 22 to 25 magarcsec -2 , BCGs are usually larger than the bulk of early 
type galaxies, and a large fraction (~ 49%) of BCGs have disky isophotal shapes. The 
differences in the scaling relations are consistent with a scenario where the dynamical 
structure and formation route of BCGs may be different from the bulk of early type 
galaxies, in particular dry (dissipationless) mergers may play a more important role 
in their formation; we highlight several possible dry merger candidates in our sample. 

Key words: galaxies : E/SOs — galaxies: cD — galaxies: cluster of galaxies — 
galaxies : photometry — galaxies 



1 INTRODUCTION 

The Brightest Cluster Galaxies (BCGs) are the most lumi- 
nous and most massive galaxies in the universe. BCGs are 
located close to the centre of the clusters of galaxies based 
on the X-ray observations or gravitational lensing observa- 
tions (e.g., Jones & Forman 1984; Smith et al. 2005). It was 
noted very early on that some BCGs show an excess of light 
('envelopes') over the de Vaucouleurs (r 1 / 4 ) profile at large 
radii (Matthews et al. 1964; Oemler 1973, 1976; Schombert 
1986, 1987, 1988; Graham et al. 1996), and a large fraction 
of BCGs are termed as cD galaxies (e.g., Patel et al. 2006). 
Although the origin of such extended envelopes is still not 
completely clear (e.g., Patel et al. 2006), the extended stellar 
halos of BCGs to surface brightness fi{r) < 25 mag arcsec -2 
are likely from BCGs themselves: the intra-cluster light has 
much lower surface brightness and only dominates at large 
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radius (r > 80kpc; for detailed discussions, see Zibetti et al. 
2005; Bernardi et al. 2007; Lauer et al. 2007). 

There has been debate on the formation mechanisms 
for BCGs and cD galaxies. Several mechanisms have been 
proposed, such as galactic cannibalism (the merging or cap- 
ture of cluster satellites due to dynamical friction, Ostriker 
& Tremaine 1975; White 1976; Ostriker & Hausman 1977), 
tidal stripping from cluster galaxies (Gallagher & Ostriker 
1972; Richstone 1975, 1976; Merritt 1985), and star forma- 
tion on BCGs by cooling flows (e.g., Fabian 1994). Therefore, 
it is important to study the photometric properties of BCGs 
using large and homogeneous BCG samples, which has now 
become available due to large surveys such as the Sloan Dig- 
ital Sky Survey (SDSS; for earlier studies, e.g. Postman & 
Lauer 1995, see references below). Moreover, the statistical 
properties of the extended envelopes of BCGs may also help 
to understand how the BCGs form and evolve, especially 
when combined with high-redshift cluster samples such as 
EdisCS (White et al. 2005). 
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Recently, there has been much progress in this area from 
both observational and theoretical fronts. Observations show 
that BCGs have distinct properties from the bulk of ellipti- 
cal galaxies. The size- luminosity relation and Faber- Jackson 
relation (L — a relation) for BCGs are significantly steeper 
than those of non-BCG elliptical galaxies (Lauer et al. 2007; 
Bernardi et al. 2007; Desroches et al. 2007; von der Linden 
et al. 2007) . It implies that BCGs are larger and have lower 
velocity dispersions than non-BCG elliptical galaxies, which 
confirms the pioneering work based on a small BCG sample 
by Oegerle & Hoessel (1991). The differences between BCGs 
and non-BCG elliptical galaxies demonstrate that the BCGs 
may form in a qualitatively different way from non-BCG el- 
liptical galaxies (De Lucia & Blaizot 2007). Furthermore, 
the steeper L — a relation for BCGs than those of non-BCG 
elliptical galaxies leads to a contradiction in the prediction 
of black hole (BH) masses using the Mbh — a or Mbh — L 
relationships. The usual Mbh — cr relation cannot predict 
black hole masses larger than 3 x 10 9 Mq in the universe, 
which is inconsistent with values estimated by other means 
(Lauer et al. 2007). 

Recent simulations and semi-analytic works in the cold 
dark matter hierarchical structure formation framework pro- 
vide a plausible picture for the formation of BCGs: they 
tend to form at high-density peaks when their inhabited 
dark matter halos collapse at high redshift, and then X-ray 
driven cooling flows allow a rapid collapse and the forma- 
tion of a stellar component (Cowie & Binney 1977; Fabian 
& Nulsen 1977; Fabian 1994). While more than half of the 
stellar mass in the BCGs may have formed before redshift 
three, BCGs still grow substantially through dry (dissipa- 
tionless) mergers when their host massive halos accrete and 
merge with other halos since redshift one (Gao et al. 2004; 
De Lucia & Blaizot 2007). This picture is largely consistent 
with observations, which does not require cooling flows to 
provide the cold gas if BCGs form late; it also overcomes a 
potential problem that the merger rate in clusters may be 
too low due to the high velocity dispersion in dynamically re- 
laxed clusters. Furthermore, dry mergers have been directly 
observed in cluster environments as well as in the field (e.g., 
Lauer 1988; van Dokkum 2005; Tran et al. 2005; Bell et al. 
2006). However, there are still no systematic investigations 
on the formation of cD-like BCGs galaxies, especially how 
the extended stellar halos of cD galaxies form: whether they 
appear during the dry merger phase, and whether they are 
the primary cause of the distinct properties of BCGs. There 
is also no agreement about the contributions of the intra- 
cluster light (e.g., Gonzalez et al. 2005; Zebetti et al. 2005; 
Lauer et al. 2007). 

The C4 cluster catalogue (Miller et al. 2005) with 748 
cluster of galaxies was constructed from the spectroscopic 
data of Second Data Release (DR2) of the Sloan Digital Sky 
Survey (SDSS). This is one of the largest published homo- 
geneous database of nearby clusters and groups (z < 0.17). 
In principle, it is suitable for statistical studies of BCGs 
(see also the more recent, much larger maxBCG catalogue 
of 13,823 clusters of Koester et al. 2007). Based on BCGs 
from this published cluster catalogue, or the un-published 
C4 cluster catalogue selected from SDSS DR3, Bernardi et 
al. (2007) and von der Linden et al. (2007) find that BCGs 
have different properties, in particular their scaling laws, 
from non-BCG early type galaxies. Their results are, how- 



ever, not entirely consistent with each other, which could 
be due to different sample selections or photometry meth- 
ods. In fact, the slope indices in the size- luminosity and 
Faber- Jackson (L — a) relations found by Bernardi et al. 
(2007) are also significantly different from those of Lauer et 
al. (2007), although the trends are the same and these two 
groups use similar methods to measure luminosities and ef- 
fective radius R e ] both are based on fitting the r 1//4 law to 
the R— band surface brightness profiles within 50 kpc or to 
some surface brightness limit of the BCGs. However, the 
samples they used are different. Lauer et al. (2007) use their 
own all sky, volume limited (z<0.05) survey of 119 BCGs 
with precise ground-based surface photometry (Postman & 
Lauer 1995). Therefore, it is worth investigating whether 
sample selections and photometry can lead to inconsisten- 
cies on scaling laws. In this paper, we construct a different, 
nearby (with redshift z < 0.1) BCG sample based on the C4 
cluster catalogue, and perform our own photometry. In par- 
ticular, we examine how the scaling relations change when 
different isophotal limits are applied, and compare the re- 
sults with those of Bernardi et al. (2007), Lauer et al. (2007) 
and von der Linden et al. (2007). 

The outline of the paper is as follows. In §2 and §3 we 
describe the selection of our local BCG sample and photo- 
metric data reduction. We present our main results in §4 
and finish with a discussion & summary in §5. Through- 
out this paper we adopt a cosmology with a matter density 
parameter Q m = 0.3, a cosmological constant Qa = 0.7 
and a Hubble constant of Ho = 70kms _1 Mpc _1 (h = 
ifo/lOOkms^Mpc -1 = 0.7). 



2 SAMPLE 

Our early type BCG sample is drawn from the early ver- 
sion of the C4 cluster catalogue (Miller et al. 2005). Given 
that BCGs are located near the centre of clusters of galaxies 
which is often crowded with galaxies and contaminated by 
intra-cluster light, to ensure the reliable photometric mea- 
surement with high enough S /N ratio for each sample BCG 
to n(r)< 25 mag/arcsec 2 that is significantly higher than 
the intra-cluster light (Zibetti et al. 2005; Lauer et al. 2007; 
Bernardi et al. 2007), we restrict our BCG sample to be 
brighter than 14.5 mag in the SDSS r-band model magni- 
tude. After we exclude 10 duplicated objects (Bernardi et 
al. 2007), there are 114 galaxies selected out of a total of 
748 BCGs in the C4 cluster catalogue. 

Due to fiber collisions in the spectroscopic data, it is 
necessary to incorporate the SDSS photometric catalogue to 
avoid missing BCGs using the C4 algorithm. The DR2 ver- 
sion of the C4 catalogue tried to correct for this by selecting 
the brightest cluster galaxy based on the photometric cata- 
log. However, as several works (e.g., Bernardi et al. 2007; von 
der Linden et al. 2007) noticed, some stars or spiral galax- 
ies (some of these are located at the edge of clusters) have 
been mis- classified as BCGs. we thus performed visual in- 
spections of two-colour images in the g-band and r-band for 
all the 114 BCG candidates and their clusters in a region 
of 200 kpc x 200 kpc and 2 Mpc x 2 Mpc, respectively. We 
found that C4 1024, C4 2092, C4 2097, C4 2031, C4 2109, 
C4 3021 and C4 3268 are seriously contaminated by adja- 
cent bright stars or foreground galaxies; C4 3047 is located 



Photometric Properties and Scaling Relations of BCGs 3 



at the edge of the corrected frame (Stoughton et al. 2002); 
C4 3235 is a star; C4 2178 and C4 3257 are nearly edge- 
on galaxies with dust lanes. In addition, there are 5 objects 
(C4 1265, C4 1332, C4 2130, C4 3098, and C4 3205) whose 
parent cluster richness, defined as the number of galaxies 
within 1 h~ x Mpc, is less than 10, and hence they may be 
wrong identifications. In the remaining 98 objects, there are 
10 BCGs that are late type galaxies with obvious spiral arms 
(C4 1023, C4 1053, C4 1324, C4 1366, C4 2001, C4 3016, 
C4 3143, C4 3246, C4 3282, and C4 3286). These are also 
excluded from our early type BCG sample, leaving us with 
88 possible early-type BCGs. We further excluded 3 objects 
(C4 1048, C4 1186 and C4 3285) with M r , 25 > -21 that are 
fainter than L* galaxies (M* ~ —21.12 for h = 0.7, see Blan- 
ton et al. 2003a) and they are unlikely BCGs (Bernardi et al. 
2007). The detailed information for all 29 rejected objects 
are listed in Table 1. Our final sample consists of 85 BCGs 
that are in the redshift mange of 0.03^z^0.09 and r— band 
(model) apparent magnitude in ~ 13.5 < mag ^ 14.5. 

Furthermore, we carefully searched for the second 
ranked cluster galaxies in the 85 C4 clusters and found that 
the BCGs are about 1.6 mag brighter on average than the 
second ranked cluster galaxies for the r— band absolute mag- 
nitude at 25 mag arcsec -2 . Only 6 BCGs have absolute mag- 
nitudes within 0.5 mag of the second ranked cluster galax- 
ies. During this process, we found one mis-identified BCG 
for C4 2020 (Abell 0119), the brightest galaxy in this clus- 
ter is UGC 00579, instead of ARK 021 as listed in the C4 
catalog. We caution that such cases may also exist in the 
unpublished C4 catalogue based on the DR3 of the SDSS 
(see von der Linden et al. 2007). 

Notice that there are 4 BCGs (C4 2049, C4 1176, C4 
3311 and C4 1035) in our BCG sample described above with 
obvious merger features. C4 2049, 1176 and 1035 are merger 
galaxies with two close nuclei and C4 3311 is with clear 
asymmetrical shapes. We will discuss these four objects in 
more detail in £14.21 

In order to perform direct comparisons with non-BCGs, 
we construct a control sample of normal early-type galax- 
ies from the SDSS DR2 spectroscopic catalogue with galaxy 
velocity dispersion measurement in the MPA gal .info cata- 
logue q These control galaxies are chosen to be in the same 
redshift and apparent magnitude ranges as sample BCGs, 
i.e., 0.03^z^0.09 and 13.5^mag^l4.5. 244 early-type galax- 
ies are selected as the control sample after excluding nearly 
edge-on objects (most of them have dust features), contam- 
inated sources, and objects identified as C4 BCGs. The dis- 
tributions of the corrected redshifts relative to the Local 
Group (see Blanton et al. 2005), (m p ,25 — m r ,2b) colour, 
and absolute magnitudes at 25 mag arcsec -2 for these two 
samples are shown in Fig. [1] It can be seen that the con- 
trol sample has a lower redshift range than those of BCGs. 
This is because non-BCGs are intrinsically fainter than the 
BCG^l (see the bottom panel), so for the same magnitude 
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2 Bernardi et al. (2007) found that all galaxies brighter than 
-24 in the r-band are BCGs. Lin & Mohr (2004) found that 
BCGs are about one magnitude (0.83) brighter on average than 
the second ranked cluster galaxies even within a small metric 
radius of 13.7/i -1 kpc. This does not contradict with the larger 
difference (1.6 mag on average) we found for our BCG sample 
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Figure 1. Distributions of corrected redshifts (top panel), 
(nig ,25 — m r ,25) colour (middle panel), and absolute magnitudes 
at 25 mag arcsec -2 (bottom panel) in the r-band for 85 sample 
BCGs (solid) and 244 control galaxies (dotted) respectively. Two 
vertical dashed lines in bottom panel show —23.5 ^ M r ^5 ^ 
—22.5, where the two samples have a significant overlap. We will 
compare scaling relations for these two subsamples in M4.4I 



limit, they reside more locally. Therefore, it is difficult to 
construct a control sample in the same luminosity range as 
that of BCGs. Instead we attempt a comparison in a narrow 
luminosity range, —23.5 ^ M r ,25 ^ —22.5 (indicated by two 
vertical dashed lines at the bottom panel of Fig. [1]) where 
there are sufficient overlaps between the control sample and 
BCGs (see £14.41 for comparisons). Furthermore, both sam- 



because we measure the magnitudes to 25 mag arcsec 2 , which 
corresponds to a larger physical size than the radius used by them. 
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pies have colour (m g ^5 — rn r ,25) ^ 0.7, implying that all of 
them are early type galaxies colour-wise (see Blanton et al. 
2003b; Shen et al. 2003 and reference therein). Also the con- 
centration index C = R90/R50 for virtually all our galaxies 
is larger than 2.6, again consistent with their being early- 
type galaxies (Shimasaku et al. 2001; Strateva et al. 2001). 
All the basic parameters for these 85 BCGs are listed in 
Table 2. 



3 DATA REDUCTION 

3.1 Estimation of the Sky Background 

Given that the SDSS image survey in the u— and z— bands 
are relatively shallow and the i— band images suffer from 
the 'red halo' effect (e.g., Michard 2002; Wu et al. 2005), we 
use the g— and r— band images in this work. The corrected 
frames have already been processed by the SDSS photomet- 
ric pipeline (PHOTO), including bias-subtraction, flat-fielding, 
cosmic ray removal and corrections for pixel defects. We use 
the astrometry for these two bands obtained by the SDSS 
astrometry pipeline (ASTR0M), which has a typical error of 
less than 0.1 (Stoughton et al. 2002). On the other hand, 
for reasons we will discuss below we perform our own pho- 
tometry. 

Although the images from the SDSS archive were pro- 
cessed with photometric corrections, as Wu et al. (2005) 
pointed out, some spurious features still exist in some im- 
ages. The spatial variation of these features was about 1-2 
ADU. Without corrections, target galaxies that happen to 
be located within such features will have inaccurate back- 
ground subtraction and poor surface brightness profiles, es- 
pecially in the outskirts of the images. We thus first correct 
for these structures following Wu et al. (2005) before con- 
structing the sky background model below. 

As Lauer et al. (2007) pointed out that the SDSS pho- 
tometric reduction systematically under-estimates the lumi- 
nosities and half-light radius of BCGs. This arises because 
the PHOTO pipeline often over-estimates the sky background 
for galaxies with large size and/or in crowded fields. Both 
problems are present for BCGs. As a result, Bernardi et al. 
(2007) performed their own photometric reductions for their 
sample. On the other hand, von der Linden et al. (2007) 
address this problem by adding up to 70% of the differ- 
ence between the local (256 x 256 pixels) and the global 
(2048 x 1498 pixels) sky background (both are available 
from PHOTO) around the BCGs to the radial surface bright- 
ness profiles and the photometry. However, the sky back- 
ground in the SDSS frames often shows spatial gradients 
and asymmetry (see the third row of Fig. [2] for an exam- 
ple). Therefore, adding a constant background value to the 
radial surface brightness profiles may not yield the most ac- 
curate photometry, at least on one-to-one basis (see Fig. 2] 
for comparisons). Therefore, we performed our own photom- 
etry following the method of Zheng et al. (1999) and Wu et 
al. (2002), which was developed to perform deep photometry 
for nearby large spiral galaxies. This method has been ap- 
plied to NGC 5907 and NGC 4565, which achieved accurate 
photometry down to an isophotal limit of ~ 29 mag arcsec -2 
in the intermediate bands of the BATC system (Fan et al. 
1996). 



To obtain accurate information of the sky back- 
ground, we first generate a 2048 x 1498 pixels (13.5' x 
9.8') background-only image using SExtractor (Bertin &; 
Amounts 1996) by masking all the detected objects with 
counts above la noise (of the whole frame) in a frame 
smoothed by a circular Gaussian with a standard deviation 
a — 3 pixels. As most of sample BCGs have sizes of a few 
square arcminutes, there are sufficient regions in the masked 
image to determine the sky background. A median filter with 
51 x 51 pixels is then convolved with the unmasked pixels, 
after which second-order Legendre polynomials are used to 
fit both rows and columns respectively (see Zheng et al. 
1999 and Wu et al. 2002). The fitted Legendre polynomials 
are then further smoothed using a circular Gaussian filter 
with a = 9 pixels to obtain our final sky background model. 
For most BCGs, we find that the sky background is tilted 
with a spatial variation about 1 — 2 ADU across the whole 
frame. We can subtract this model from the frame to obtain 
the sky-free image. Our sky subtraction procedure is simi- 
lar to that in Gonzalez et al. (2005). Furthermore, the sky 
background subtracted not only includes the contributions 
of the intra-cluster light (which is at a much lower level), but 
also from other astrophysical sources (dust emission etc., see 
Gonzalez et al. 2005 for a discussion). 

Our procedure is illustrated in Fig. [2] for C4 3086. The 
first, second and third rows of Fig. [2] show the frame, the 
masked frame and the smoothed sky background in the 
r— band. Notice that the sky background shows a gradient 
across the frame. The bottom panels show the distributions 
of counts in the sky-subtracted frame for all unmasked pixels 
and for the local vicinity around the target BCG. If the sky 
background model is successful, then we expect the back- 
ground counts to follow a Gaussian distribution with a mean 
close to zero. This is indeed the case, as can be seen from the 
bottom panels of Fig. [2] The distributions in the r— band for 
the whole image and the local region are well-fit by Gaus- 
sians with dispersion of 5.52 ADU and a mean value of about 
-2 x 10" 2 ADU. 

After background subtraction, we geometrically align 
the two-colour frames, and then trim the corrected frames 
to 501x501 pixels centred on the target. We run SExtractor 
again on the trimmed frame to generate a c SEGMENTATION ; 
image, which identifies all objects with flags in the frames. 
A mask image with all the detected objects except the 
galaxy of interest flagged can then be obtained from the 
'SEGMENTATION 5 image. We carefully examined all the mask 
images in two colours and corrected a few bad images manu- 
ally to create good mask images for all galaxies. Photometry 
is then performed on the trimmed images with the masked 
areas excluded from the reduction; we discuss the details 
below. 

3.2 Isophotal photometry 

To compare the scaling relations, such as the size-luminosity 
relation and Faber- Jackson relation, of BCGs with those of 
the bulk of elliptical galaxies, we need to perform accurate 
surface photometry. In this work, we measure the isophotal 
magnitudes within a certain radius, instead of magnitudes 
based on a model (e.g., a de Vaucouleurs [1948] model or a 
Sersic [1968] model) as a large fraction of BCGs is not well 
described by such simple models. Furthermore, to see how 
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Table 1. Reasons for the rejection of 29 out of 114 C4 BCGs. 



C4 ID 


R.A.(J2000) 


Dec.(J2000) 


reason for rejection 


(1) 


(2) 


(3) 


(4) 


1023 


1 664897 

-L kJkJ . UU4:0^ 1 




idle iy JJC 


1053 


998 ^D7Qfi1 


4 987^^^ 


Idle iy|Jc 


1324 


99*3 D4^Qfi8 


-D 9^fi1 08 


Idle iypt? 


1366 






Idle iy|Jc 


2001 


^^D 8fi9fi74 


1 4 ^9^841 


Idle lyjJc 


3016 


1 87 4^9774 


R4 n^ni ^ 


Idle iy|Jc 


3143 


947 1 D4DRQ 


41 1fi84fi^ 


Idle lypc 


3246 


1 ^fi ^7998^ 


^fl D8Q9QD 
ou.uoyzyu 


Idle ly|Jc 


3282 


94Q ^1 8497 


44 41 8997 


Idle lyjJc 


3286 


1 84 9*39^1 1 


fi^ 4DQQ99 


Idle iy|Jc 


1024 


99fi R88n^1 


-1 9^1 70Q 
- j_ . J. i uy 


A V "1 T 1 1~ 1 T t 1 1 Y 1 Q "f" O /i I \"\ r Kri ry n"i~ C 1" Q T' 

dJllldlllllldltJd L»y Ulljilll oldl 


2031 


^94 77*W\k 


D 7DR1 Q1 
-u. t uuiyi 


contaminated by bright star 


2092 


94 99^1 ^8 


_ U. JOOU 1 -L 


CUIlldillllldlcCI L)y Dllglll oldl 


2097 


^8 7Q71 fifi 

OOO. 1 a t 1UU 


1 HAQ^OQ 
- j. . uiyjzy 


contaminated by bright star 


2109 


94 ^1 4n*S1 


-Q 1 Q7fi1 1 
-a . iy 1 Ull 


CUIlldillllldlcLl L>y Dllgill oldl 


3268 


oqq qqqi p;n 


^4 R71 1 D^i 


contaminated by bright star 


3021 


247.489689 


40.630707 


contaminated by foreground spiral 


2178 


351.297310 


15.199853 


edge-on galaxy with dust lane 


3257 


249.074493 


44.135700 


edge-on galaxy with dust lane 


3235 


119.610850 


37.732302 


a star 


3047 


158.245415 


56.748148 


at the edge of the frame 


1265 


168.680626 


4.024536 


richness < 10 


1332 


167.820450 


-0.824176 


richness < 10 


2130 


50.701678 


-6.694635 


richness < 10 


3098 


246.907102 


42.638233 


richness < 10 


3205 


135.367430 


55.044549 


richness < 10 


1048 


147.895859 


1.112059 


L < L* 


1186 


183.990581 


3.305910 


L < L* 


3285 


259.895020 


56.630066 


L < L* 



the scaling relations vary as a function of isophotal limits, 
we measure the photometric parameters to four different 
surface brightness limits, 22, 23, 24, and 25 mag/arcsec 2 . 
We emphasize that the photometry for the control sample 
has been performed in the same manner as that for BCGs. 

The surface photometry analysis is performed following 
Wu et al. (2005). We briefly outline the procedures below; 
the readers are referred to that paper for more details. The 
IS0PH0TE/ELLIPSE task in IRAF is used to fit each of the 
trimmed background-subtracted images with a series of el- 
liptical annuli from the centre to the outskirts. The width of 
annuli is chosen to increase uniformly in logarithmic steps, 
with the semi-major axis radius increasing by 10% between 
two adjacent annuli. The size of the annuli in the outer parts 
is therefore larger, which suppresses the shot noise in the 
outer regions where the signal-to- noise ratio (S/N) is lower. 
The ellipticity, position angle, and other quantities are also 
fitted simultaneously. The r— band images are first used to 
define the isophotal annuli, which are then applied to the 
g— band images. The coordinates of the photometric peak in 
both r— band and g— band are obtained by the DAOPHOT 
package and are fixed during fitting. 

We estimate the seeing of our galaxies in two colours 
by averaging the star profiles in the corresponding frames 
respectively. The average values of the FWHMs of the see- 
ings in both colours are ~ 1.6", which are similar to the 
measurements by the SDSS collaboration. We integrate the 
observed surface brightness profiles directly to estimate the 



isophotal apparent magnitudes and half-light radius Rso, 
within which half of the integrated flux of galaxies is con- 
tained, under each isophotal limit. The absolute magnitude 
(M) is derived from the obtained isophotal apparent mag- 
nitude (m) by M = m — 5 log(.Dz,/10 pc) — A — k, taking 
into account the extinction (A) in each filter by the SDSS, 
and the k-correction (using the KC0RRECT algorithm of Blan- 
ton & Roweis 2007). Other photometric parameters, such as 
ellipticities, coefficients of lowest order deviations from per- 
fect ellipse ci4 and 64 (see Bender et al. 1988), are derived 
following Hao et al. (2006). Notice that the a^/a parame- 
ters are weighted by the intensity between twice the see- 
ing radius and the radius where the surface brightness is 
25 mag arcsec -2 in the r— band. We discuss the isophotal 
shapes in Sect. l4~3l 

All our radial profiles in this paper use the equivalent 
radius of an ellipse, y/ab, where a and b are the lengths of 
semi-major and semi-minor axes of the ellipse. The surface 
brightness profiles in the two bands are performed account- 
ing for the Galactic extinction correction, cosmological dim- 
ming, and the k-correction (Blanton & Roweis 2007). The 
photometric errors are estimated as in Wu et al. (2005). 
The observational errors in the surface brightness profiles in 
each band include random errors (e.g., readout noise, the 
shot noise of the sky background and BCGs) and the error 
from the sky subtraction. 

To illustrate the performance of our photometry, we 
show the comparison of our determined Petrosian parame- 
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ters with those of SDSS in Fig. [3] One can see that the SDSS 
pipeline under-estimates the sizes (and thus luminosity) of 
BCGs, which has been pointed out by other workers (e.g., 
Lauer et al. 2007). There is a systematic trend: the under- 
estimation becomes more serious for larger (brighter) BCGs. 
This is in good agreement with the results of Desroches et 
al. (2007, see their Fig. 1.). 

As mentioned before, von der Linden et al. (2007) 



adopted a simple prescription to correct for the sky sub- 
traction problem. In Fig. [4] we show a comparison between 
their surface brightness profile, ours and that from the SDSS 
for three example galaxies. It can be seen that our results 
are in very good agreement with those of von der Linden et 
al. (2007). While there are small differences at faint surface 
brightness (25 mag arcsec -2 ), such differences have little im- 
pact on their results since they adopted a surface brightness 
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of 23 mag arcsec -2 to measure the physical quantities. So 
for statistical purposes, their method can offer an efficient 
way of performing photometry, particularly at high surface 
brightness limits. 



4 RESULTS 

4.1 The Extended Envelopes of BCGs 

As discussed in the introduction, the surface brightness pro- 
files of the majority of cD galaxies show strong deviation 
from a perfect de Vaucouleurs or Sersic profile due to cD 
galaxies embedded in an extensive luminous stellar halo. In 
order to quantitatively measure the extended envelope of 
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Figure 2. An example of sky subtraction (for target BCG C4 
3086) in the r— band. The first row shows the SDSS corrected 
frame (2048 x 1498 pixels; 13.5' x 9.8'), the second row shows the 
masked frame, and the third row shows the fitted sky frame. The 
inner box in each frame has 501 x 501 pixels (3/3 x 3/3), and the 
outer box has 751 x 751 pixels (4/7 x 4/7). The bottom panels 
show the normalised distributions of flux in the unmasked region 
after the sky background subtraction, for the whole frame (left) 
and the local region between the two boxes (right) in the r— band. 
The solid lines in the two panels are the best-fit Gaussians to the 
distributions. 




R Petro50 [ k P C ] 

Figure 3. Comparison between our measured Petrosian half- 
light radii to those provided by the SDSS for our BCG sample. 
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Figure 4. Comparison between the surface brightness profiles 
obtained using our method (thick solid lines) and those from the 
SDSS (dotted lines) and von der Linden et al. (2007, thin solid 
lines) for three example galaxies in the r— band. The vertical and 
horizontal dashed lines indicate twice seeing radius and one per- 
cent of the sky brightness respectively. Notice for the left and 
middle galaxies, our results almost overlap with those from von 
der Linden et al. (2007). 



BCGs, we use an objective method based on the Petrosian 
rj(r) profiles (Patel et al. 2006), defined as 



r](r) = Li(r) - (/x(r)}, 



(1) 



where ji(r) is the surface brightness in magnitudes at radius 
r and is the mean surface brightness within r (Pet- 

rosian 1976). There is a distinct signature of a plateau in 
the Petrosian 77 (r) profiles for cD galaxies with an extended 
stellar halo (see below). Such a plateau is not present for 
normal elliptical galaxies which are usually well-fit by a de 
Vaucouleurs' surface brightness profile (e.g., Kjaergaard et 
al. 1993; Brough et al. 2005). Furthermore, notice that the 
plateau is not present for the Sersic surface brightness pro- 
file that better fits some ellipticals (see Fig. la of Graham 
et al. 1996). 

Fig. [5] shows the surface brightness profiles and Pet- 
rosian j]{r) profiles for 3 BCGs. For BCG C4 1187, its profile 
is well fit by a Sersic with n « 4.53 (in the g-band) model, 
and the 77 (r) profile does not show any distinct feature. In 
contrast, for C4 3034 and C4 3077, their surface brightness 
profiles deviate significantly from a single Sersic profile (or 
the r 1//n profile, n — 4 for the de Vaucouleur profile), and 
the 77 (r) profile shows a plateau. Furthermore, we can see 
that the more significant the deviation of the surface bright- 
ness profile from the a single Sersic profile (i.e., the more 
prominent the stellar halo in the outskirt), the flatter the 
plateau in the r\ (r) profile becomes. 
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Figure 5. Images, surface brightness, r](r) and dr](r)/d log(r) 
profiles for three typical examples BCGs (C4 1187, C4 3034 and 
C4 3077). The images shown in the top panels are in the r— band 
and have a size of 501 x 501 pixels. In each image, north is up 
and east is to the left. The second row shows the radial surface 
brightness profiles in the g— band (lower curve) and r— band (up- 
per curve) respectively. The solid line is the best-fit Sersic model 
for each surface brightness profile that is obtained by fitting inten- 
sity profile from twice seeing radius to radius at 25 mag arcsec -2 
respectively. The vertical and horizontal lines are twice the see- 
ing radius and one percent of the sky brightness in the g— band 
(dashed lines) and r— band (dotted lines) respectively. The third 
row shows the observed (dots) and smoothed (line) Petrosian 77(7*) 
profiles in the r— band. The bottom row shows the correspond- 
ing dr](r)/d log(r) profile. The vertical solid line in each panel 
indicates the minimum position, 



To quantitatively determine how a surface brightness 
profile deviates from the r 1 ^ 71 profile we calculate the gra- 
dient of rj(r), j(r) = drj(r)/d\og(r). As mentioned before, 
a single r 1 ^ model profile always results in a monotonic 
7(7-) profile for any n. The bottom panels of Fig. [5] show 
this quantity as a function of radius. It is clear that galaxies 
with the most significant deviation from the model profile 
also have deeper valleys in the 7(7*) profile. We label the 
radius of the minimum in 7(7*) as R m i n , and also define a 
depth, Dmin, as the difference between the minimum and 
maximum of 7(7*) outside twice the seeing radius and i? m in- 

The right histogram in Fig. [6] shows the distribution 
of D m in for our sample BCGs (the median value is ~ 0.8). 
We can see that the distribution of D m in shows a weak bi- 
modality, but is still consistent with a continuous distribu- 
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Figure 6. Correlation of D m i n with i^mm, r— band absolute 
magnitude, and size for 85 BCGs, measured out to an isophotal 
limit 25 mag arcsec -2 . The right-most panel shows the histogram 
of the D m in distribution together with the Poisson errors; the 
median value is ~ 0.8 (indicated by a thick dotted line). 



tion within Poisson error bars. This shows that the devia- 
tion of the surface brightness profiles from the a single Sersic 
profile does not have any sharp transitions. Therefore, it is 
difficult to give an unambiguous criterion to separate cD 
galaxies from non-cD BCGs. 

The left panel of Fig. [6] shows the correlations between 
D m in and R m \ n . The null- hypothesis that there is no correla- 
tion between D m i n vs. R m in can be rejected at the 7.5 x 10 -4 
level. So statistically the larger the radius i? m in, the more 
significant the extended envelope a BCG appears to have. 
The middle and right panels in Fig. [6] show the correlations 
between D m in with the r— band luminosity and the half-light 
radius for BCGs. The null-hypothesis that there are no cor- 
relations between D m in vs. M r ,25, and D m in vs. Rso,25 can 
be rejected at the 3.8 x 10 -5 , and 9.6 x 10 -6 confidence lev- 
els. These correlations demonstrate that more luminous and 
larger BCGs tend to have more extended stellar envelopes 
and therefore the fraction of cD galaxies increases. 

4.2 Mergers and Envelopes in BCGs 

As mentioned in 32 that there are four galaxies in our BCG 
sample that have either two close nuclei or other clear merg- 
ing signatures, such as broad faint fans or asymmetry in the 
morphology. The top panel of Fig. [7] shows the image of 
C4 1176, from which one can clearly see two nuclei and a 
broad surrounding fan to the north of the image. The SDSS 
provides the spectra for both nuclei, which are shown in the 
bottom panels of Fig. [7] From the spectra the redshifts of the 
two nuclei are 0.0738 and 0.0736 respectively, corresponding 
to only 60 kms -1 difference in the line-of-sight velocity. The 
angular distance between the two nuclei, ~ 4.7", implies a 
projected separation of just ~ 6.7 kpc. Therefore, C4 1176 is 
most likely a merger galaxy with two nuclei. Furthermore, 
the spectroscopic properties of these two nuclei are similar; 
both are typical of early type galaxies with weak emission 
lines of LINER characteristics, indicating that BCG C4 1176 
may be a dissipationless (dry) merger system. 

For the other three BCGs, we can see from the top 
panel of Fig. [8] that C4 2049 and C4 1035 also have two 
possible close nuclei. From the SDSS spectra, we find that 
the projected distance between the two nuclei for C4 1035 is 
only ~ 1.7 kpc and the relative line-of-sight velocity is 120 
kms -1 , and thus may be another merging BCG. For C4 
2049, the SDSS has only spectrum for one nucleus. The an- 
gular distance between these two nuclei is ~ 5.5". If they 
are at the same redshift, then the projected distance is 
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Figure 7. The image of merging BCG C4 1176 with 
two nuclei shown in the top panel. North is up and east 
is to the left. The bottom two panels show the spectra 
of the two nuclei (the line symbols are from SDSS, see 
http://www.sdss.org/gallery/gal_spectra.html); the redshifts of 
the two nuclei are almost the same, indicating they are merg- 
ing. Notice the broad low-surface brightness fan to the north of 
the image. 



~ 7.7 kpc. Hence C4 2049 is another possible merger BCG. 
For C4 3311, although one cannot see double nuclei at its 
centre (down to the SDSS spatial resolution), it shows clear 
asymmetry in shape, which is indicative of merging activi- 
ties. 

The second, third and bottom rows of Fig. [8] show the 
surface brightness profiles, rj{r) profiles and dr/(r) / d\og(r) 
profiles for all four merger/possible merger BCGs. As can 
be seen, there are obvious plateaus in the 77 (r) profiles for 
3 out of 4 merger BCGs (C4 2049, 1176, and 3311). The 
Dmi n values measured from the drj(r) /d\og(r) profile are 
1.34, 1.34, and 1.04, for C4 2049, 1176, and 3311 respectively, 
which are much larger than the median D m in value (~ 0.8) 
for the whole BCG sample. Therefore, these three BCGs 
clearly have extended envelopes in their outskirts and are 
likely cD galaxies. 

The high fraction of cD galaxies in merger BCGs and 
the fans corresponding to cD galaxy envelopes suggest that 
the extended stellar halos of BCGs are from mergers and 
they appear shortly after merging takes place. 

In fact, the role of dry (dissipationless) mergers in the 
formation of ellipticals has recently received much atten- 
tion. As pointed out by Bell et al. (2006) and van Dokkum 
(2005), the typical morphological signatures of dry mergers 
are broad stellar fans, short tidal tail with similar colour 
as the galaxy itself, and asymmetries at very faint surface 
brightness levels. The four examples of merger BCG candi- 
dates in our sample shown at Fig. [8] indeed have all these 
signatures, but they also show even more direct merger signs: 
in three of these four cases, there are two close nuclei. Fig. [9] 
shows ten more examples of possible merger systems from 
our sample. They have broad low-surface brightness fans or 
distinct asymmetries (top panels); many visually appear to 
inhabit in environment with multiple small satellites (bot- 
tom panels). These images provide strong evidence for the 
connection of mergers with the formation of stellar halos 
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Figure 8. Images, surface brightness, 77(7*) and dr](r)/d\og(r) 
profiles in the r— band for four likely merger BCGs in our sample. 
The format is the same as in Fig. [5] In each image, north is up 
and east is to the left. 




Figure 9. The panels in the top row show the images of BCGs 
with obvious dry merger signatures, such as broad fans or faint 
asymmetry shapes. In the bottom row, all 5 images of BCGs are 
surrounded by many possible satellites or are in pair or triple 
systems; these are possible signatures of interaction and mergers. 
The values of D m in for all these dry merger candidates are larger 
than the median value (0.8) in the Z} m in distribution for sample 
BCGs. In each image, north is up and east is to the left. 



in BCGs. The 59 BCGs with spectroscopy from SDSS have 
spectra similar to the examples shown in Fig. [7] indicat- 
ing our BCG samples are typical early-type galaxies. Their 
spectra do not show significant star formation at the centre, 
although one third of our BCGs show weak emission lines 
with LINER characteristics. Therefore, the mergers appear 
to be largely dissipationless (dry). We will return to the role 
of mergers in the formation of BCGs in a future study. 
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Figure 10. Structure parameter a^/a versus r— band absolute 
isophotal magnitude, M r ,25, for the BCGs (solid circles) and the 
control sample (open circles) respectively (left panel). The median 
error bar in 0,4/0,, and the fractions of galaxies with boxy isophotes 
for BCGs and the control sample are shown in the inset. The right 
panel shows the histograms of the 04/0 distributions for the BCGs 
(solid histogram) and the control sample (dotted histogram). The 
median values are shown as solid and dashed lines respectively. 



4.3 Isophotal Shapes 

We also investigated the shapes of isophotes for BCGs and 
compared with those of the control sample. The isophotes of 
ellipticals are usually well fit by ellipses. However, small but 
significant deviations exist, which are usually described by 
the Fourier components of the deviations. The most signifi- 
cant one is the a^/a (Lauer et al. 1985; Bender et al. 1989), 
with a positive (negative) cu/a indicating a disky (boxy) 
isophotal shape. 

Fig. 1 101 shows how the cla/cl parameter varies with the 
luminosity of galaxies for both the BCGs and the control 
sample. It is clear from Fig. [10] that a^/a decreases as the 
luminosity of galaxies increases, hence the fraction of boxy 
galaxies increases for both samples; this trend has been 
pointed out by many authors (e.g., Bender et al. 1989; Hao 
et al. 2006). However, the boxy fraction for BCGs (51%) is 
not much larger than that of the bulk of ellipticals (37%), 
even though BCGs are more luminous than the bulk of el- 
lipticals in the control sample. In particular, we can see 
from Table 2 and Fig. \10\ that the boxy fractions of BCGs 
with luminosities in the range of —24 < mag ^ —23 and 
mag< —24 are ~ 55% and ~ 65% respectively. So there 
are still a large fraction of very luminous early type galaxies 
with disky isophotes (cu/a > 0), which is consistent with 
the HST image results on BCGs by Laine et al. (2003). 
Faber et al. (1997) and Lauer et al. (2005) pointed out 
that galaxies brighter than My — —22 mag (correspond- 
ing to M r ~ —22.6 mag, after correcting the difference in 
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Figure 11. Luminosity (L r ), size (R50), mean surface brightness 
within R50 ((/i5o)) and the concentration factor (C' = R90/R50) 
as a function of the isophotal limit. The solid and dotted lines 
denote the BCG and control samples respectively. Error bars are 
the 68.3% confidence levels of the median value. 



H , and taking V - r = 0.3 ± 0.1, Krick et al. 2006) tend to 
have 'cored' luminosity profile, while galaxies fainter than 
My — —20.5 mag (corresponding to M r « —21.1 mag) tend 
to have 'power-law' luminosity profiles (see also Rest et al. 
2001; Lauer et al. 2005). Other works (e.g., Rest et al. 2001) 
show that bright galaxies are slowly rotating, boxy systems 
(Davies et al. 1983) while faint elliptical galaxies are fast 
rotating, disky systems (Bender 1988; Bender et al. 1989; 
Nieto et al. 1988). However, BCGs appear to somewhat con- 
tradict with this trend: some very luminous BCGs still have 
disky isophotes. The flattening in elliptical galaxies can arise 
either from an-isotropic motions or from rotation. Faint el- 
lipticals appear to be flattened by rotation (e.g., van den 
Bosch et al. 1994; Lauer et al. 2005), while the luminous 
ones may be supported by an-isotropic motions. As we have 
shown in the last subsection, BCGs may have experienced 
fairly recent merging events, and thus may still retain some 
of the orbital angular momentum from merging, and their 
isophotes may be disky. 



4.4 Scaling relations 

There have already been many works on the scaling relation 
of BCGs based on homogeneous BCG samples (Bernardi et 
al. 2007; Desroches et al. 2007; von der Linden et al. 2007; 
Lauer et al. 2007). However, so far no work has investigated 
how the slopes change with the isophotal limit, which can in 
principle provide more knowledge of the dynamical structure 
of BCGs. 



4-4-1 The size-luminosity relation 

The left panels of Fig. 1111 show how the size (R50) and lu- 
minosity (L r ) vary with the isophotal limits. It is clear that 
the size (R50) and luminosity of BCGs increase more sharply 
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Figure 12. The size-luminosity relation for BCGs (red solid cir- 
cles) and the control sample (blue open circles) for four isophotal 
limits, 22, 23, 24, and 25 magarcsec -2 respectively. The median 
error bar for the data points is indicated in the bottom right cor- 
ner of each panel. The black solid and dashed lines are the best 
power-law fits for the BCGs and control sample respectively. The 
power-law index ex (i?5oocL a ) is shown at the top of each panel. 

than that of the control sample as the isophotal limit be- 
comes deeper. The right panel shows how the mean surface 
brightness within i? 50 ((mso)) and the concentration factor 
(C' — R90/R50) vary with the isophotal limit. We can see 
that (/iso) of BCGs decreases more sharply than that of 
control sample, but C = R90/R50 of BCGs has the same 
trend as that of the control sample as the isophotal limit 
becomes deeper. Therefore, BCGs appear to have signifi- 
cantly larger average sizes and are more diffuse in surface 
brightness. These trends agree with those found in previous 
studies (e.g., Lugger 1984; von der Linden et al. 2007). 

Fig.ll2lshows the size-luminosity relations for our BCGs 
and the control sample, measured to four isophotal limits, 
22, 23, 24 and 25 mag/arcsec 2 respectively. The solid and 
dotted lines in each panel show the best power-law (R50 oc 
L a ) fits; the a value is indicated at the top left corner of 
each panel0 

Clearly the size-luminosity relation becomes steeper as 
the isophotal limit becomes fainter for both the BCG and 
the control sample. However, we can also see that the size- 
luminosity relation for BCGs is steeper than that of the 
control sample at the same isophotal limit. In particular, the 



3 Notice that the intrinsic scatters in the data points are much 
larger than the statistical error bars, and so the formal % 2 per 
degree of freedom is large, indicating the presence of systematics. 
To approximately account for these, we renormalised the statisti- 
cal errors by a constant factor until the % 2 per degree of freedom 
is unity. The errors on the power-law indices are obtained using 
these renormalised errors. A similar procedure was adopted in 
von der Linden et al. (2007). 



size-luminosity relations for BCGs is already slightly steeper 
(but not very statistically significant) than those of the bulk 
of ellipticals when the isophotal limit is 22 mag arcsec -2 , cor- 
responding to a physical half-light radius of about 5 kpc (see 
Fig. [TT]) . It indicates that the deviation of size- luminosity re- 
lation for BCGs from that of bulk of ellipticals may already 
occur for the inner regions of ellipticals, and thus is not be- 
cause most BCGs have extended stellar halos. Rather, this 
suggests the deviations arise from different dynamical struc- 
tures of BCGs from the bulk of ellipticals. Bernardi et al. 
(2007) arrived at a similar conclusion by comparing the size- 
luminosity relations for BCGs with and without extended 
stellar envelopes (see also von der Linden 2007). 

The power-law index (0.88±0.04) for the size-luminosity 
relation we derived for BCGs at an isophotal limits of 25 
mag/arcsec 2 is slightly smaller than the value (0.92) ob- 
tained by Bernardi et al. (2007), but the power- law index 
(0.74 d= 0.03) for the control sample at the same isophotal 
limit is a bit larger than their value (0.62). The difference 
in the slope measured by us and by Bernardi et al. (2007) 
may be because we measure the isophotal magnitudes by 
integrating the observed surface brightness profiles directly, 
while Bernardi et al. (2007) measured the magnitude using 
a de Vaucouleurs or Sersic model for the surface brightness 
profile. In fact, when we use the same method as Bernardi 
et al. (2007) to perform photometry, the power-law index 
for the size-luminosity relation we obtain is almost the same 
(0.93 ± 0.03) as that of Bernardi et al. (2007), indicating 
that different ways of photometry can give different power- 
law indices in the size-luminosity relation. 

von der Linden (2007) also derived the size-luminosity 
relation down to 23 magarcsec -2 for their BCGs. The slope 
they derived is a = 0.65 =b 0.02, similar to that for their 
control sample (a = 0.63 =b 0.02). Our slope for the BCGs 
(0.77=b0.04) is steeper than their value at the same isophotal 
limit, while the slope for our control sample (a — 0.65±0.03) 
agrees well with theirs. We have 66 overlapping galaxies with 
the BCG samples used by von der Linden et al. (2007). If we 
use only these galaxies, we derive a slope similar to theirs. 
Our steeper slope appears to arise from the fact that our 
BCGs are brighter since their sample includes more fainter 
BCGs which makes their slope shallower (von der Linden 
2007, private communication). Therefore, the sample con- 
struction also influences the power-law index in the size- 
luminosity relation. 

The power- law index a — 1.18 d= 0.06 in the size- 
luminosity relation by Lauer et al. (2007) for their core 
galaxies with My ^ —21 (most core galaxies with My ^ 
—22 are BCGs) measured in a similar way as Bernardi et 
al. (2007) is significantly steeper than those of Bernardi et 
al. (2007), von der Linden (2007) and our results. It again 
illustrates that how samples are constructed and the meth- 
ods used to do the photometry lead to different power-law 
indices in the size-luminosity relation. Comparisons between 
different studies should be aware of this. 

4-4-2 The Faber- Jackson relation and Fundamental plane 

It is well known that the luminosities and velocity disper- 
sions of many early type galaxies satisfy the Faber- Jackson 
relation, L oc a ^ , where (3 ~ 4 (Faber & Jackson 1976). 
Given that the SDSS spectra were taken within a 3" fiber 
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Figure 13. The Faber- Jackson relation for the 59/85 BCGs and 
244 control sample galaxies with spectroscopic information for 
four different isophotal limits (22, 23, 24, and 25 mag arcsec -2 ) 
respectively. The symbols are the same as in Fig. 1121 The power- 
law indices (3 are shown at the top of each panel for the BCGs 
and the control sample. The median error bar in the velocity 
dispersion is shown at the bottom right in each panel. 

aperture (J0rgensen et al. 1995; Bernardi et al. 2003a), we 
correct the velocity dispersions for BCGs and the control 
sample following von der Linden et al. (2007, see their equa- 
tion 4). The uncorrected value of a for each BCG is shown 
in Column 5 of Table 2. 

Fig. [13] shows the L — a relation for 59 out of 85 
BCGs and the control sample with SDSS spectral informa- 
tion, measured to four different isophotal limits, 22, 23, 24, 
and 25 magarcsec -2 . Clearly at a given surface brightness 
limit, the slope for BCGs in the Fabor- Jackson relation is 
steeper than that of the bulk of early type galaxies. The 
index (3 for both the BCG and control samples increases 
as the isophotal limits become deeper. The /3 value for the 
control sample is around 4, consistent with the canonical 
value for the Faber- Jackson relation. However, this is not 
the case for BCGs: only when the isophotal limit is 22 
magarcsec -2 , corresponding to the inner region of BCGs, 
f3 is about 4. For deeper isophotal limits, the f3 values for 
BCGs are significantly larger, approaching 5.6 for a surface 
brightness limit of 25 mag arcsec -2 . von der Linden (2007) 
found f3 = 5.32 ± 0.37 for their BCGs at an isophotal 
limit of 23 mag arcsec -2 , which is within 2a of our value 
((3 = 4.52 ± 0.41). For their control samples, they found 
(3 = 3.93 ± 0.21, within la of our result [fl = 3.98 ± 0.19). 
Lauer found /3 = 6.5 =b 1.3 for their core galaxies with 
My < — 21, which is significantly larger than those of values 
mentioned above. Note that the sample Lauer et al. (2007) 
used is not a pure BCG sample, although most of their core 
galaxies with My < — 22 are BCGs. 

The steeper power-law indices in the size-luminosity 
and (to a less extent) in the Faber- Jackson relations may 
arise because BCGs are brighter than the control sample 
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Figure 14. The size-luminosity relation (left) and Faber-Jackson 
relation (right) for the BCG and control sample galaxies at 
—23.5 ^ M rj 25 ^ —22.5. The top panels use all 94 control galax- 
ies and 45 BCGs in this luminosity range. The bottom panels are 
the results for one realisation for 45 BCGs and randomly selected 
45 control galaxies with the same luminosity distribution. The 
symbols are the same as in Fig. [12] The power-law indices are 
shown at the top of each panel for each sample. 



or because there are intrinsic differences between BCGs 
and the bulk of elliptical galaxies. In order to differenti- 
ate these two possibilities, we repeat the above analysis for 
our BCG and control samples in the same magnitude range, 
-23.5 < M r ,25 < -22.5 (see 22). Fig. [14] shows the size- 
luminosity (left) and Faber-Jackson (right) relations. The 
top panels use all the 94 control-sample galaxies and 45 
BCGs within the luminosity range. Since Fig. [1] shows the 
control sample is still fainter on average than the BCG sam- 
ple within this narrow range, we have randomly selected 45 
out of the 94 control-sample galaxies so they match well 
with the luminosity distribution of the 45 BCGs. The sta- 
tistical results for different realisations are similar; the bot- 
tom panels shows the results for one realisation. From this 
figure, we clearly see that the power-law indices for the size- 
luminosity relations for BCGs are steeper than those for 
bright non-BCG elliptical galaxies; the power-law indices 
for the Faber-Jackson relations are also steeper, but only at 
~ la levels, and so are not statistically significant. To sum- 
marise, the data are consistent with a scenario where the 
cluster environment plays an important role on the BCG 
formation. However, the numbers of galaxies in our control 
and BCG samples are still too small to differentiate the two 
possibilities mentioned above. 

Fig. 1151 shows the fundamental plane for our BCGs and 
the control sample. While they appear to have different 
power-law slopes in the Faber-Jackson relation, both BCGs 
and the control sample follow roughly the same fundamen- 
tal plane, although the scatters in the BCGs appear smaller. 
At bright isophotal limits, ji — 22, and 23 mag arcsec -2 , the 
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Figure 15. Fundamental planes for the 59 BCGs with spectra 
(dots) and 244 control sample (open circles) for four different 
isophotal limits respectively. The solid line in each panel is the 
fundamental plane determined by Bernardi et al. (2003b) using a 
uniform sample of ~9000 early type galaxies. The median error 
bars are shown at the top left. 

fundamental planes for BCGs and the control samples ap- 
pear to differ somewhat. These results are in good agreement 
with Bernardi et al. (2007). 



5 SUMMARY 

In this paper we have studied the properties of 85 BCGs 
carefully selected from the C4 cluster catalog (Miller et al. 

2005) . To address the problems of the SDSS pipeline PHOTO 
in crowded fields and for large galaxies, we have performed 
our own photometry for these BCGs, as well as for a control 
sample of galaxies in similar redshift and apparent magni- 
tude ranges. Particular attention was paid to accurate ex- 
traction of the sky background in clusters using a method 
developed for high-precision photometry for large, nearby 
spiral galaxies (Zheng et al. 1999; Wu et al. 2002, see g)). It 
came as somewhat a surprise to us that the simple recipe of 
von der Linden et al. (2007) gives surface brightness profiles 
similar to ours, and hence for many statistical purposes their 
simple recipe can be adopted as a more efficient alternative. 

We analysed the surface brightness properties of our 
BCGs using the r](r) profiles (see eq. [I]). We demonstrate 
that its gradient, 7(7*) = drj(r) / d\og(r) , can be used as an 
indicator of extended envelopes of cD galaxies (Patel et al. 

2006) . Furthermore, we find that the depth of valley, D m i n 
in the 7(7*) profiles can be used to quantitatively measure 
the properties of the plateaus in the surface brightness pro- 
files. Generally, the deeper the valleys are (corresponding to 
a larger value of .Dmin), the natter and wider the plateaus 
are in the 77 (r) profiles. However, the D m [ n parameter ap- 
pears to vary in a continuous way with the galaxy param- 
eters, such as their size and luminosity. While there is a 



clear trend that more luminous and larger BCGs have larger 
stellar halos. the continuous variation makes it difficult to 
un- ambiguously classify cD or non-cD galaxies based on the 
surface brightness profiles alone. 

We measured the photometric properties of galax- 
ies to four different isophotal limits, 22, 23, 24, and 
25 mag arcsec -2 , and used these to investigate the scaling 
laws and their dependence on these limits. The latter is 
of particular relevance as some previous investigations de- 
rived scaling relations for BCGs using different samples and 
different isophotal limits (e.g., the study of von der Lin- 
den et al. 2007 used an isophotal limit of 23 mag arcsec -2 ). 
We find that the size- luminosity relation and Fabor- Jackson 
(L — a) relations of BCGs are consistent with being steeper 
than those of the bulk of early type galaxies for a control 
sample selected within the same absolute magnitude range 
(—23.5 < M r ,25 < —22.5), although the latter is not statis- 
tically significant due to the limited sample sizes. 

Furthermore, as the photometric limit becomes deeper, 
the size-luminosity relation for BCGs and the bulk of 
the ellipticals becomes steeper. However, there are already 
(small) differences between BCGs and regular ellipticals at 
22 mag arcsec -2 , a surface brightness that probes the inner 
parts of galaxies. This, together with the relatively small 
difference in the scaling relations between 24 mag arcsec -2 
and 25 mag arcsec -2 , suggests that there may be intrinsic 
differences in the dynamic structures between BCGs and 
the bulk of early type galaxies (Bernardi et al. 2007; Lauer 
et al. 2007). One possibility may be that the amount of 
dark matter in the central parts of BCGs and regular ellip- 
ticals is different (von der Linden et al. 2007), or perhaps 
the relative importance of an-isotropy and rotation may be 
different. In this context, it is interesting to note that we 
find a large fraction of very luminous early type galaxies 
with disky isophotes. This appears to somewhat contradict 
with the usual expectation that the most luminous galaxies 
should show boxy isophotes. Also notice that it is difficult to 
measure high-redshift galaxies to very faint surface bright- 
ness due to the cosmological dimming. Thus a proper com- 
parison between local and high-redshift samples will need to 
explicitly account for the surface brightness dependence we 
found in the paper. 

In this work we uncovered four strong dry merger candi- 
date BCGs (see Fig. [8}, three of which have two close nuclei 
within lOkpc and with nearly identical radial velocities. We 
also found another 10 BCGs in our sample with possible dry 
merger signatures (see Fig. [9}, as advocated by Bell et al. 
(2006). This strongly suggests merging at the centre of clus- 
ters is an important physical process that affects a variety of 
properties of BCGs, both in photometry and kinematics. An 
analysis of the statistics for the fraction of dry mergers in 
BCGs is ongoing and the results will be reported elsewhere. 
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